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Abstract 

We present a temperature dependent photoluminescence study of silicon optical nanocavities 
formed by introducing point defects into two-dimensional photonic crystals. In addition to the 
prominent TO phonon assisted transition from crystalline silicon at ^ 1.10 eV we observe a broad 
defect band luminescence from ^ 1.05 — 1.09 eV. Spatially resolved spectroscopy demonstrates 
that this defect band is present only in the region where air-holes have been etched during the 
fabrication process. Detectable emission from the cavity mode persists up to room-temperature, 
in strong contrast the background emission vanishes for T > 150 K. An Ahrrenius type analysis 
of the temperature dependence of the luminescence signal recorded either in-resonance with the 
cavity mode, or weakly detuned, suggests that the higher temperature stability may arise from an 
enhanced internal quantum efficiency due to the Purcell-effect. 
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Due to its indirect electronic bandgap, crystalline silicon is rarely used as the active 
emitter in semiconductor optics [Ij. Interband light emission is predominantly a phonon- 
assisted process and silicon has, therefore, a very poor internal quantum efficiency. However, 
the development of efficient silicon based light emitters would pave the way toward CMOS 
compatible monolithic optical interconnects and, therefore, signal processing speeds much 
higher than currently provided by silicon micro-electronics [2 [3l HJ [5l [6] . 

Recently, enhancement of the photoluminescence (PL) intensity from crystalline silicon at 
room-temperature has been observed using two-dimensional (2D) silicon photonic crystals 
(PhCs) with photonic point defect nanocavities [Zll^. Such defect PhC nanocavities modify 
the spatial emission profile of light; due to the in-plane photonic band gap a much larger 
fraction of luminescence is emitted perpendicular to the slab than parallel to it [9j . In addi- 
tion, the high Q-factors that are attainable using PhC nanocavities [121[II1[I2] together with 
their small mode volumes may lead to an enhancement of the internal quantum efficiency 
of the active material due to an increase of the radiative emission rate via the Purcell-effect 
[13]. Other work demonstrates enhanced PL emission from internal light emitters embedded 
in photonic crystal slabs [HI [15]. 

In this letter, we present a detailed investigation of the spectrum and temperature sta- 
bility of the PL emission from crystalline silicon PhC nanocavities. Analysis of our results 
indicates that both the phonon satellites of the interband silicon emission and surface defect 
states are responsible for the luminescence of PhC cavity modes. Most interestingly, we 
can detect luminescence from the cavity mode up to room temperature whilst the back- 
ground emission intensity rapidly reduces below our detection threshold for T > 150 K. An 
Arrhenius type analysis of the temperature dependent data, either in resonance with the 
cavity mode or spectrally detuned from it, indicates that the local radiative emission rate is 
enhanced in resonance via the Purcell-effect. 

The samples investigated were fabricated from silicon-on-insulator (SOI) wafers with a 
d = 250 nm thick active silicon layer on top of a 3 |Lim thick layer of Si02. A 2D PhC is 
patterned into the upper silicon layer as illustrated in fig. [l](a). We defined a triangular 
lattice of air holes with a period of a = 275 nm in the silicon slab as shown by the scanning 
electron microscope (SEM) image in the left panel of fig. |l| (b). This was done using: electron- 
beam-lithography and subsequent SF6/C4F8 reactive-ion-etching. These techniques allow us 
to control the radius r of the air holes with a precision of ±2 nm. As a final fabrication 
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FIG. 1: (a) Left panel: schematic cross-sectional representation of the photonic crystal nanocavity 
structures investigated. Right panel: layer sequence in the active region, (b) Left panel: SEM 
image showing a photonic crystal from top. Right panel: zoom-in to the modified L3 defect where 
the outer holes are shifted by 0.15 lattice constants. 

step the underlying Si02 is removed by hydrofluoric acid. Low mode-volume nanocavities 
{VMode ^ 1 • i^/^T) were realized by omitting three air-holes in a row and by shifting the 
lateral holes away from the cavity center by 0.15 lattice constants to form modifled L3 cav- 
ities [TBJ. A typical scanning electron micrograph (SEM) image of the samples investigated 
is shown in the right panel of flg. [l] (b). 

Spatially-resolved optical measurements were performed using a micro-photoluminescence 
(//PL) spectroscopy setup. The sample was excited by a diode-pumped CW solid-state-laser 
emitting at Xhaser = 532 nm. We focussed the laser beam using a 50x microscope objective 
(NA = 0.5) to a spot size of 1 |am. The resulting PL signal was collected through the 
same objective and dispersed by a 0.32 m imaging monochromator equipped with a 600 
lines/mm grating and a liquid nitrogen-cooled InGaAs linear diode array. 
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In fig. [2] (a) we present room-temperature //PL spectra recorded from a series of L3 
nanocavities as the normalized radius r/a is increased from 0.28 (bottom) to 0.34 (top). 
The excitation power density used in these measurements was 550 kWcm~^. Using these 
measurement conditions the intensity of the emission from the unpatterned region of the 
device was below our detection sensitivity. In strong contrast, each of the spectra recorded 
from the PhC nanocavities clearly reveals five distinct peaks (marked with arrows on fig. [2] 
(a)) that shift to higher energy with increasing r/a ratio [19j. The extracted peak energies 
are plotted in fig. [2] (b), clearly demonstrating that all peaks (labeled Ml to M6) shift in a 
similar way to higher energies with increasing r/a, as previously reported in refs. 0|8]. In 
fig. [2] (c) we compare the measured PL spectrum for r/a = 0.28 with the spectral positions of 
mode-emission obtained from photonic bandstructure simulations (vertical red lines). Using 
the geometric parameters extracted from SEM images, we obtain six eigenmodes from the 
simulation results. The insets in fig. [2] (c) show the calculated electric field distributions that 
reveal three modes with longitudinally oriented field distributions (M1,M2,M6) and three 
modes with vertically oriented field distributions (M3,M4,M5). In full accordance with 
Iwamoto et al [7J, we observe the strongest peak intensities for modes that have a high 
vertical extraction efficiency (M1,M3,M6). Simulation and experiment for the longitudinal 
modes are generally in very good agreement. A discrepancy is observed between simulation 
and experiment for M5, which is not observed in the spectra, probably due to a higher 
sensitivity of mode Q and frequency to structural disorder. This tends to lead to a low-Q 
and, thus, broad emission from M5. The comparison between simulation and experiment 
clearly shows that we observe PhC mode-emission in the spectral vicinity of the silicon 
interband luminescence and that we are able to clearly distinguish the cavity mode emission 
from the background. 

In fig. [3] (a) we compare low temperature (T = 22 K) /xPL spectra recorded close to 
the Ml mode emission from three different spatial positions as indicated on the inset. The 
measurement positions were: on the nano-cavity (blue trace), on the PhC away from the 
cavity (red trace) and on the unpatterned material (green trace). The peak at ^ 1.10 eV 
arises from the TO-phonon replica of the interband emission from the silicon [17j . In addition 
to the TO-line, we observe a broad emission band when detecting on the PhC that extends 
from ^ 1.05 — 1.09 eV on the low energy side of the TO-replica . In bulk silicon one would 
not expect emission in this range at low temperatures, since all phonon replica are sharply 
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FIG. 2: (a) Room-temperature PL spectra recorded from a series of L3 PhC cavities with different 
normalized air-hole radii ranging from r/a — 0.28 (bottom) to r/a — 0.34 (top), (b) Detailed 
analysis of the five peaks (labeled Ml to M6, marked with arrows in (a)) shifting to higher energy 
with increasing r/a. (c) Comparison between experimental PL data (black curve) and simulated 
spectral mode positions (vertical red lines) of the six L3 modes for r/a — 0.28. Insets show 
calculated electric field profiles of the eigenmodes. 

defined [17j. We investigated the intensity of the broad emission band as a function of the 
excitation position on the sample. Selected results of these measurements are presented in 
fig- [3] (b) where we plot the integrated PL intensity from LOG eV to LOS eV (highlighted 
orange in (a)) for a spatial scan across the PhC (compare with fig. [3] (a)-inset). Here, we 
clearly observe a direct correlation between PL intensity of the emission in the LOG — LOS eV 
band and the number of air-holes in the probed region. This strongly indicates that the 
broad emission band arises from phonon-mediated recombination from surface defects on 
the side walls of the air holes created by the RIE etching process. A large number of surface 



5 



Energy (eV) Position (jjm) 

1.04 1.06 1.08 1.10 1.12 4 8 12 16 20 

-r— I — ' — I — ' — I — ' — r 




CO 

"cD ^ 
O 

CO 



•(b) 

" / \ 


1 1 1 1 1 

\ - 




II 1 II III! hil li iljiili^^^^^^^^^^ 



CQ 

— ? 

I— K 


Q- 



1000 



500 



I I I I I I I I I 



1 ' 1 ' 1 ' 1 ' 1 ' 1 ' 1- 

T = 24 Ki 



reference PL: 
r/a = 0.25 




— ? 
p- 

c 

C/) 



Tl 
|— 

5" 

I— K 





o 
■a 

C/) 



1.05 1.06 1.07 1.08 1.09 1.10 1.11 1.12 

Energy (eV) 

FIG. 3: (a) PL at T = 22 K for three different detection positions (indicated by the inset) on 
the sample structure, (b) Integrated PL of broad emission band at the low-energy side of the TO- 
replica (range highlighted orange in (a)) when scanning across the PhC. Inset shows correlation with 
etched air-hole structure, (c) Red data: integrated mode intensity of M3 as a function of emission 
energy for r/a ranging from 0.25 to 0.20. Black line: reference PL spectrum for r/a = 0.25 at 
T = 24 K. 



defect states with different trapping energies would lead to the inhomogeneously broadened 
emission band on the low energy side of the TO-phonon replica. To investigate the origin of 
the emission in the PhC cavity modes we controUably shifted M3 through the TO-phonon line 
and the surface defect band by progressively decreasing the air-hole radius from r/a = 0.25 
to r/a = 0.20 and extracted the intensity of the M3 emission. The results are summarized 
in fig. |3] (c), where we plot the integrated intensity of M3 as a function of peak position (red 
data points). For comparison, we plot the PL spectrum recorded from the r/a = 0.25 sample 
(black line). The mode intensity clearly tracks the spectrum of the TO replica and follows 
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also the intensity trend of the surface defect band. These observations indicate that the 
cavity modes are pumped predominantly via the phonon replica and also more weakly via 
the surface defects introduced by the fabrication process. This arises from the fact that the 
luminescent defects are close to the air-silicon interface, where the electric field amplitude 
of the eigenmodes as shown in fig. [2] (c) are weak, leading to a weak coupling between mode 
and surface defects. 

We continue to analyze the temperature stability of the PL spectrum recorded from the 
PhC cavity site. Typical temperature dependent data recorded from the Ml mode emission 
for a PhC with r/a = 0.34 from 24 K to room-temperature (295 K) is presented in fig. [3] (a). 
We observe a clear and rapid reduction in the intensity of the TO-phonon replica, the Ml 
mode emission and the surface defect band, respectively. This is caused by an increasing 
importance of non-radiative decay channels, like Auger and free carrier recombination [18j. 
For further quantitative analysis we describe the temperature dependent PL intensity Ipl{T) 
as a function of the external photon detection probability P, the excitation rate R and the 
internal quantum efficiency r]int{T)'. 

IpL{T) = P-R-r,,nt{T), (la) 

with r]int{T) = y . . (lb) 

J- Rad + -1- NR[-L ) 
-0 _/^~^^ 



and Tnr{T) = P^^ • exp l^J^ j ^ (Ic) 
In these equations Tjiad is the interband radiative recombination rate and Tnr^T) is the 
non-radiative recombination rate. The temperature-dependence of T^RiT) is described by 
distributing excitation amongst the two levels, |2) and |1) separated by Ea according to 
Boltzmann statistics. A schematic few level-diagram of the system is presented in fig. [3] (b). 
State |1) is pumped at an excitation rate R and decays either via a radiative channel to |0) or 



non-radiatively over |2). Eqn. ^a^to eqn. [Tclare directly obtained from rate equations gov- 



erning the steady state populations of |1) and |2) using the approximation P^^ ^ P^^(T). 
This condition ensures that the non-radiative channel dominates carrier recombination at 



elevated temperature as expected for silicon. By reformulating eqn. la we obtain 



where Ipl{T^^) is the extrapolated PL intensity for T ^ K [20j. In fig. [i] (c) we plot 
\ii{IpL{T^^) / Ipl{T) — 1), as a function of inverse temperature in an Arrhenius-type graph. 
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FIG. 4: (a) PL spectra as a function of temperature for a PhC with r /a — 0.34. (b) Schematic 
level diagram illustrating the decay from state |1) via a radiative channel to state |0) and via a 
non-radiative channel over state |2). (c) Arrhenius-type analysis for the TO-replica (black squares), 
defect band emission (blue circles) and Ml mode emission (red triangles). The activation energies 
Ea extracted from linear fits label the three curves 
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according to eqn. [2| The analysis was performed by spectrally integrating over the TO- 
replica (black squares), the surface defect band PL from 1.07 meV to 1.08 meV (blue circles) 
and the Ml cavity mode emission (red triangles) [21]. For all datasets we observe a straight 
line as predicted by eqn. [2] supporting the vailidity of our analysis. From the slopes of these 
curves we extract an activation energy Ea of 8 ± 2 meV (surface defects), 10 ± 2 meV (TO- 
replica) and 10 ± 2 meV (Ml mode) for non-radiative recombination processes. To obtain 
these fits we only took into account Ipl{T) for T < 100 K where sufficiently intense signal 
from mode, defect band and TO-replica are available to make a valid comparison. From 
the similar values of Ea for the cavity mode and TO-replica we conclude that the mode 
is predominantly excited via the TO-replica, as already deduced from the data shown in 
fig. 1(c). 

Comparing mode and TO-replica in fig. [4](c), we see from eqn. [2] that the lower values 
of \ii{IpL{T^^) / Ipl{T) — 1) from the mode emission can only be explained via the term 
^Rad/^%fR) since both emission peaks exhibit the same activation energy. Hence, the data 
suggests that the presence of the cavity mode results in a ratio of Tpiad/^%R which is larger 
than for the TO-replica. As shown in refs. [7, 9j, an enhanced PL signal can also be caused 
by a redistribution of the spatial emission profile via the photonic crystal structure. However, 
this effect does not increase the radiative recombination rate. Therefore, the results indicate 
that the enhanced photoluminescence from the cavity mode may be caused by an enhanced 
internal quantum efficiency r]int{T) due to a larger radiative carrier recombination rate caused 
by the Purcell-effect [22]. This would explain the observed temperature stability of the mode 
emission up to room-temperature, in strong contrast to the vanishing background emission 
for T > 150 K. 

In conclusion, we presented a temperature dependent investigation of PL in Si PhC 
nanocavities. We suggest two mechanisms being responsible for the luminescence of cavity 
modes, namely phonon-mediated recombination from charge carriers in the electronic band 
states and recombination from charge carriers trapped in surface defect states. Furthermore, 
we observed an enhanced internal quantum efficiency in spectral resonance with the cavity 
mode emission. 
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